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We investigate wind tunnel turbulence generated by both conventional and multi-scale 
grids. Measurements were made in a tunnel which has a large test-section, so that possible 
side wall effects are very small and the length assures that the turbulence has time to 
settle down to a homogeneous shear-free state. The conventional and multi-scale grids 
were all designed to produce turbulence with the same integral scale, so that a direct 
comparison could be made between the different flows. Our primary finding is that the 
behavior of the turbulence behind our multi-scale grids is virtually identical to that 
behind the equivalent conventional grid. In particular, all flows exhibit a power-law decay 
of energy, ~ t~^, where n is very close to the classical Saffman exponent of n = 6/5. 
Moreover, all spectra exhibit classical Kolmogorov scaling, with the spectra collapsing on 
the integral scales at small k, and on the Kolmogorov micro-scales at large k. Our results 
are at odds with some other experiments performed on similar multi-scale grids, where 
significantly higher energy decay exponents and turbulence levels have been reported. 



1. Introduction 

In recent years a number of wind-tunnel experiments and numerical simulations have 
focused on quasi-homogeneous turbulence generated by multi-scale grids. (See, for ex- 



ample, Hurst & Vassilicos 2007 and Nagata et al. 2008 ) When the grids concerned 



possess certain fractal-like properties, this is sometimes referred to as fractal-generated 
turbulence. In some of these experiments the turbulence appears to behave in unex- 
pected ways. For example. Hurst & Vassilicos (2007) report unusually high values of 
Re\ — u\/u, perhaps twice that associated with the equivalent conventional grid. (Here 

A is the Taylor micro scale, — (u^) / (^{du^ / dx)^^ , v the viscosity, and u = ^ (u^) /3 

the characteristic fluctuating velocity.) The same paper also reports energy decay char- 
acteristics, such as power-law decay, (u^) ^ , with very high decay exponents (n ~ 2 
rather than the classical ?! ~ 1.2 — > 1.4), and also exponential decay rate for some grid 
geometries. This is surprising as there are good arguments why classical fully-developed 
homogeneous turbulence can decay no faster than n ~ 10/7 (see e.g. Appendix [A] of this 
paper, |Krogstad fc Davidson) |2010[ or [Ossai fc Lesieur[ |2000[ for a summary of these 
arguments). The implication is that turbulence generated by multi-scale excitation can 
exhibit non-classical characteristics, which has led to a re-examination of the assumtions 



underlying energy decay laws (see for example George & Wang 2009 and Hosokawa 



2008). 



However, this suggestion is somewhat at odds with the evidence of direct numerical 
simulations in periodic cubes, where it is usually found that, after some transient, the 
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behavior of the turbulence is largely independent of the precise form of the initial energy 
spectrum (see, for example, [Ossai fc Lesieur , ^2000| ) . Indeed, it is often argued that, once 
fully developed, all the turbulence remembers of its initial conditions is the prefactor, Cm, 
in the expression E{k — )■ 0) = Cmfc"*, Cm being an invariant for to < 4 (Ossai & Lesieur 



2000 Davidson 2004 Ishida et al. 2006 and Davidson 2009). 



There is, however, a different interpretation of the experiments of [Hurst fc: Vassilicos| 
(2007), one which is consistent with the evidence of the numerical simulations. Immedi- 
ately behind a multi-scale grid the turbulence is far from homogeneous, and there is a 
significant production of turbulent energy due to transverse gradients in the mean veloc- 
ity. Moreover, this inhomogeneity in the mean and fluctuating velocities is stronger than 
for a conventional grid and it persists for many tens of integral scales downstream of the 



mesh (Nagata et al.. 2008). This might not matter if all measurements are made well 



downstream of a multi-scale grid, where approximate homogeneity is recovered. In con- 
ventional grid turbulence it is considered good practice to restrict measurements to, say, 
a; > 25 — >■ 30M in order to avoid the inhomogeneous region behind the grid (e.g. Comte- 



Bellot & Corrsin 1966 Gad-el-Hak & Corrsin 1974). It is possible, therefore, that the 



apparently non-classical results obtained for some of the grid geometries of [Hurst fc Vas 



silicos ( 2007 ) are a consequence of a distinct lack of homogeneity immediately behind the 



grid, i.e. a transient phenomenon which disappears further downstream. Indeed, [Hurst fc] 
[Vassilicos 2007[ themselves note that more extensive tests are required in order to assess 
the role of inhomogeneities close to the grid. 

In order to test the hypothesis that turbulence behind multiscale grids will develop 
in a classical way (after an initial transient), we have conducted a series of wind tunnel 
experiments using two multi-scale grids of the cross class described in [Hurst fc Vassilicos| 
(2007), and compared these results with data obtained using a conventional grid. In 
order to ensure that the comparison is meaningful, the dimensions of the grids were 
chosen such that the integral scale of the turbulence some short distance downstream of 
each grid, £oi is virtually the same in all three cases (to within a percent or so). Crucially, 
the experiments were carried out in a large working section, and measurements taken 
over the extended range 80^o < a; < 400^0i which translates to 50 A/ < a: < 240M for 
the conventional grid. Excellent homogeneity was obtained throughout the measurement 
range for all three grids and the long working section allowed us to determine the energy 
decay characteristics with considerable accuracy. As we shall see, the fully-developed 
turbulence generated by all three grids is virtually identical. In particular, there is no 
significant difference in the behavior of Re\, and in all cases the decay exponent n in 
the expression (u^) 



t " is very close to the classical prediction of 



Saffman 



n = 6/5. (For a discussion of decay exponents in classical turbulence see e.g 



2004 and Krogstad fc Davidson 2010 



( [1967[), i.e. 
Davidson 



2. The Experimental Set-up 

The experiments were performed in the large recirculating wind tunnel described in 



Krogstad fc Davidson (2010). The tunnel test-section has transverse dimensions of 2.7m 



x 1.8m (measured at the start of the test section) and is 12m long. There is an adjustable 
roof to compensate for the growth of the sidewall boundary layers and the grids were 
mounted upstream in the test section contraction to improve isotropy. From the location 
of the grid to the entrance of the test section, the area contraction ratio was 1.48 and 
the test section starts x — 1.2m downstream of the grid. 

All three grids were produced from 2mm thick sheet metal. The conventional grid 
(labeled eg) has square holes 30mni x 30mm punched at 40mm spacing, giving a mesh 
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Figure 1. The two multi-scale grids used in this study. 



size of AI = 40mm, a bar width oi t — 10mm, and a solidity oi a — 44%. The tests on 
this grid were all performed at a Reynolds number of ReM = UM /v = 3.6 x 10'', where 
U = 13.5m/s was the mean speed in the tunnel. 

The first of the multi-scale grids (labeled msgl) is similar to the cross-grid-type(a) of 



Hurst & Vassilicos ( 2007 ) (which we shall label eg — a) and is shown in Figure 1 (a) It 



has bar widths ranging from ti — 8mm down to ^3 — 2mm, and mesh sizes ranging from 
Ml = 6Amm to M3 = 15mm. The solidity of msgl is also a = 44%. These measurements 
were taken at U = 14.0ni/s. 

The second multi-scale grid {msg2) is shown in Figure [1(b) As for msgl, the bar 
widths vary from ti — 8mm down to — 2mm, though the mesh sizes are larger, with 
All — 88mm to AI3 = 21mm. This reduces the solidity of msg2 to a — 33%. This grid 
was tested a.t U = 15.5m/s. 

As we shall see, the turbulence produced by these grids becomes more-or-less homo- 
geneous and fully developed at around x = 2m, at which point the Kolmogorov mi- 
croscale is 77 0.22mm — s- 0.26mm. On the other hand, the integral scales at a; = 2m, 
io = i{x = 2m), obtained in the usual way by integrating the longitudinal correlation 
function 

{u^{x)u^{x + r)) ^ 

{uD ' ^'-^^ 

turn out to be ^Q = 23.9mm (for eg), £0 = 23.6mm (for msgl), and £q — 23.4mm 
(for msg2), respectively. (The measurements of 77 and i are discussed in f|5]) Note the 
uniformity of £q across the grids. Note also that the geometric length-scales associated 
with the two multi-scale grids almost span the range of dynamic scales associated with 
the turbulence, from around 9r] up to several integral scales. Finally we note that, in 
terms of £0, the tunnel cross-section is approximately 115^o x 80£o, thus ensuring that 
there is minimal influence of the side-wall boundary layers. Measurements are taken up 
to a considerable distance downstream of the grid, allowing an accurate estimate of the 
energy decay exponent, n. 

The data was obtained using single and two component hot-wire anemometry. For 
the measurements of the decay of ^(ua,)^^, spectra, the length scales 77 and £ etc. a 
purpose made 2.5/.tm partly etched Platinum-10%Rhodium straight single wire probe 
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was produced. The active wire length was 0.5mm. For two component measurements 
and for global checks such as spanwise homogeneity etc. an X-wire probe with ±45° 
nominal wire angles was produced. For this probe bfim partly etched wires were used 
with 1mm wire lengths. The two wires were also separated in the spanwise direction by 
Az = 1mm. By careful handling of the two delicate probes it was possible to maintain 
the same wires throughout the whole experiment. Hence there should be no effects in the 
data that can be attributed to changes in probe characteristics. 

The probes were operated at an overheat temperature of about 320 degrees using in- 
house manufactured anemometers which were tuned to a frequency response of at 
least 20kHz. {f^ is here calculated from = (1.3t)~^, where r is the time for the 
anemometer output to recover to within 10% of the asymptotic level after a square wave 
step change has been inserted at the top of the Whetstone bridge in the anemometer.) 
The output from the anemometers were suitably offset and amplified to span as much as 
possible of the ±10 volt range of the 16 bit acquisition card used before the signals were 
passed through an AC coupled Krohn-Hite amplifier and low-pass filter unit. The filter 
frequency fc was set by inspecting the dissipation spectra of a few initial measurements 
with very high filter settings, fc was then set at the frequency where noise first started to 
affect the dissipation spectra and a new set of data was obtained at a sampling frequency 
which was slightly higher than 2/^. 

For verification purpose only, a number of two component laser Doppler anemometry 
(LDA) measurements were taken for each grid, but these data were not used in the data 
analysis presented here, since LDA tends to produce too high turbulence levels as the tur- 
bulence intensity drops below about 1%. Further information about the instrumentation 
and data analysis may be found in Krogstad & Davidson (2010). 



3. Tests for Homogeneity and Isotropy. 

Let us start by documenting the relative degree of homogeneity and isotropy in the 
experiment. Figure [2] shows the time-averaged centre-line velocity, normalized by Uq = 
U{x = 2m), for all three grids as a function of x. For x > 2m, U is constant to within 
±0.5%, with Re = Uio/u = 2.1 x lO'', 2.2 x 10* and 2.4 x 10'', respectively, for the three 
grids. However, for x < 2m there is a noticeable streamwise acceleration due to the effect 
of the contraction extending slightly into the first part of the test section. In the light 
of the streamwise acceleration, we shall ignore data for x < 2m (i.e. x < 80£o)- Figure [2] 

1/2 

also shows the streamwise variation of Re\ = (m^) X/i/ for all three grids. The present 
values are only about one third of those reported by Hurst & Vassilicos ( 2007 1 and there 
is no significant difference between the conventional and multi-scale grids. 

1/2 

The spanwise distributions of the turbulence intensity, T„ — u/U, where u = (u^) , 
measured at a; « 2m and x 6m are shown in Figure |3j along with the spanwise 
distribution of the skewness, Su, and fiatness, F^^ of Bi x ^ 2m. The measurements 
were made over a span corresponding to about z ss ±4M, where M represents the largest 
mesh of each grid. Since M is only a well defined scaling length for the conventional grid, 
we have chosen to scale z with £q, hence the apparent differences in the measurement 
ranges in the figure. 

Again, all three grids are seen to behave in a similar manner, with variations in 
limited to about ±3% at x = 2m and ±2% at x = 6m. The skewness and flatness 
distributions are also well behaved in all three cases, with Su being slightly negative and 
Fu ~ 2.95, values which are typical of most grid turbulence experiments. We conclude 
that, for X > 2m, the turbulence is relatively homogeneous. 
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Figure 2. Streamwise distributions of Rex (left axis) and the test section centre-line speed, U, 
(right axis). U is normalized by C/ at a; = 2m. 



The streamwise development of the anisotropy of the turbulence, as measured by the 
ratios {ul)/{ul), {ul) / {uf) and (^) /■^{ul), where (g^) = (u^) ^ {ul+ul + ul), 
was measured for each grid (Figure B| . For all cases it was found to follow closely the 
distrib utions reported for the conventional grid in Figure 3 of [Krogstad fc Davidson] 
(2010). We show here just the three streamwise developments of (g^) /3(w^) which is 
seen to be very similar for the 3 cases and the ratio is very close to the isotropic value of 
1. For all grids (u^) / (w^), (u^) / (ul) and (^q^) /3 (u^) were close to unity 'di x — 2m, 
but (w^) / (uy) and (u^) / (u^) were found to divert slowly with increasing x. The largest 
departure from isotropy was observed at the exit of the test section, where (u^) / (m^) 
reached values in the range of 0.8 to 0.9. There was a corresponding growth in (m^) / (uy), 
to produce the almost constant (q^) /3 (m^) ratios shown in the figure. 

Figure [5] shows the skewness and flatness distributions of the streamwise fluctuations 
measured along the centre line of the test section. The flatness is constant throughout 
the measurement range at Fu ~ 2.95 for all grids. The skewness is slightly different for 
the grids initially, as was shown in Figure [3'(b)[ but they all tend steadily towards S'„ — ^ 
as the flow develops downstream. 

In most of what follows we shall restrict the analysis of our data to the region x > 2m, 
to avoid any acceleration effects caused by the inhomogeneities in U, and to a; < 8m, 
because of the increased levels of noise for large x as the turbulence level drops well 
below 1%. For completeness, however, the full set of data will be included in the plots so 
that the reader may see where and how the data departs from classical decay behavior. 
From the data presented above we see that, in this range, the turbulence is reasonably 
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Figure 3. Spanwise distributions of; (a) at a:: ~ 2m (left axis) and x ~ 6m (right axis), and 
(b) skewness, Su, of at (left axis) and flatness, Fu, of Ux at x — 2m (right axis). 



homogeneous and that, although there is some anisotropy, the levels are not excessive 
and are comparable to, if not better, than in most other experiments. 



4. The Energy Decay Rate. 

We shall now estimate the decay exponent, n, in the decay law (u^) ~ for all 
three grids. First, however, perhaps it is worth saying something about the values this 



exponent might be expected to take on theoretical grounds. As noted in Davidson (2004 



2009) and Krogstad & Davidson (2010), there are two classical predictions for n. One 
arises when E{k — >■ 0) ~ fc^, a situation called Saffman turbulence, and the other when 
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Figure 4. Streamwise distributions of (g^) /3(ii^). 

E{k 0) ^ fc**, so-called Batchelor turbulence. In the former case it may be shown that 
L = J (u • u') = constant , (for E{k) ~ k'^) , (4.1) 



where (u • u') is the usual two point velocity correlation (see Saffman 1967), while in 
the latter case 



/ = - y (u • u') dr = constant , (for E{k) ^ k^) , 



(4.2) 



with the caveat that Eq. ( |4.2[ ) holds only once the turbulence is fully- developed (see Ishida 
et al. . 2006 ) . These integrals are dominated by the large scales and so self similarity of 



the large scales (when it applies) demands 

u'^i^ — constant , (for Saffman turbulence) , 

u'^i'^ — constant , (for Batchelor turbulence) . 
These may be combined with the empirical, but well supported, law 



dt ' T 



A — constant 



(4.3) 
(4.4) 

(4.5) 



to give n — 6/5 (Saffman's exponent; see Saffman 1967 ) in E{k) ■ 
n = 10/7 (Kolmogorov's exponent; see Kolmogorov 1941 ) in E{k) 



fc^ turbulence, and 
k'^ turbulence. The 



first of these exponents was observed in the experiments of Krogstad & Davidson (2010) 



and the second in the numerical simulations of Ishida et al.\ (2006). While other decay 
exponents have been proposed from time to time, it is natural to keep these two classical 
predictions in mind when examining the experimental data. However, there is a slight 
complication which arises when comparing these predictions with experiments: in wind 



0.6 



0.4 



S 0.2 







-0.2 







P.- A. Krogstad and P. A. Davidson 



# A 



. ^ H " ^ „ 



X (Conventional grid) 
A (Multi scale grid #1 ) 
V F (Multi scale grid #2) 



; (Conventional grid) 
□ (Multi scale grid #1) 
o S (Multi scale grid #2) 



OO o 



o 



OOo 



oOo 



o oO 



8 



x(m) 



2.8 



2.6 F„ 



^2.4 



2.2 



10 



Figure 5. Streamwise distributions of (g^) /3(ii^). 



tunnel data the coefficient A can vary slowly along the test section, and this causes slight 
departures from the ideal values of n = 6/5 or n = 10/7, even when Eq. (4.3) or (4.4) 



hold true (see Krogstad & Davidson 2010). Indeed, we shall see shortly that just such a 
slow variation of A occurs in our experiments. 

Let us now estimate the decay exponents for our three grids by comparing the exper- 
imental data with the power law 



C/2 



X — Xq 
4 



(4.6) 



where, as before, is the integral scale &t x — 2m. Perhaps we should note from the 
outset that it is notoriously difficult to obtain reliable estimates of n. There are several 



reasons for this. First, as shown in Krogstad & Davidson (2010), the virtual origin xq 
in Eq. (4.6) does not correspond to the point where the turbulence ffi'st matures, but 



rather is located somewhat upstream of that position, and we do not know in advance 
where xq will lie. Second, if the range of x/i^ is too short, the decay exponent becomes 
very sensitive to the choice of the unknown xq. Third, if data from the inhomogeneous 
region close to the grid is included in the fit, then higher (and misleading) values of 
n are usually obtained. Fourth, for large x/Iq the turbulence intensity is low, so that 
noise starts to become problematic, obvious sources being electronic noise and pressure 
fluctuations arising from the side-wall boundary layers. In summary, then, to obtain 
reliable estimates of n we require that: (i) xq be chosen systematically and with care; (ii) 
there is an extended x/-^o range; (iii) anomalous data close to the grid must be excluded; 
and (iv) data for large x/£o should be ignored when the noise level becomes excessive. 

Figure^shows (u^) /t/^ and (g^) /St/^, obtained using two different alignments of the 
x-wires, all plotted as a function of x. This log-log plot demonstrates that, for all three 
grids, there is a clear power-law relationship between and x. It also shows that (u^) 
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Figure 6. Streamwise distributions of (it^) /i7^ (filled symbols) and (g^) /3(7^ (open symbols). 



and (g^) follow essentially the same power law. Instead of attempting a direct power-law 
fit to this data, it is convenient to rewrite Eq. (4.6 1 as 



In 



t/2 



= Ina — nln 



Xo 



4 



(4.7) 



to which a linear fit may be made. In view of the difficulty in obtaining reliable estimates 
of n, three different fitting procedures will be used and their results compared. All three 



methods are described in detail in Krogstad & Davidson (2010). 



4.1. Regression method 
Here we search for the fit that gives the smallest variance, cr^, between the data and Eq. 



(4.7 



(4.7 



. The procedure works as follows. Starting with all of the data, a linear fit to Eq. 
is made and noted. Data at the most downstream location is then removed one 
by one until a minimum in is obtained. This then determines Xmax, the point beyond 
which noise starts to have a significant effect on the data. With x^ax fixed, a similar 
procedure is then applied at the most upstream location and used to determine Xmin , the 
location of the start of the admissible data. The best fit values of a, xq and n for these 
values of Xmin and Xmax are then taken to be the optimum values for these parameters. 
The exact locations of Xmin and Xmax vary slightly from grid to grid and between fitting 
to (u^) and (9^), but typically Xmin ~ 1.5m (i.e. ~ 60^o) and Xmax ~ 8w (i-S- ~ 330£o)- 
This procedure applied to (m^) yields n = 1.13, 1.11, and 1.25 for grids eg, msgl and 
msg2 respectively, as shown in Figure [7] The corresponding virtual origins are at xq = 
0.26m, 0.43m and 0.30m. Marginally different exponents were found when applying this 
method to (g^), and the range of values obtained using regression are noted in Table 1, 
alongside the range of exponents estimated using the two procedures described in the 
following sections. 




regression 
method 



regression 
method 



local exponent 
method 



n 

majdmum 
decay range 



1 Conventional grid | 


0.26m 


1 1.13±0.02 1 


1.14±0.02 


1 1.17±0.04 1 


1 Multiscale grid 1 | 


0.43m 


1 1.12±0.02 1 


1.17±0.02 


1 1.19±0.03 1 


1 Multiscale grid 2 | 


0.30m 


1 1.25±0.02 1 


1.25±0.02 


1 1.23±0.03 1 



Table 1. Estimates of the energy decay exponent n obtained using three diflterent methods. Note 
that the uncertainty in n indicated above for a given method does not represent the statistical 
uncertainty in n, but rather the range of values obtained using a given procedure for a randomly 
chosen set of data. 



4.2. Local exponent method 

Measurements of {u^} were made at intervals of Ax = 200mm or Aa; w 8.5£o with a 
positioning accuracy better than ex = 1mm and so local exponents can be estimated 



using two-point differencing applied to Eq. (4.7 1 



In 






In 


(x-\-Ax—xq) 




{x — Ax—xq) 



(4.8) 
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Figure 8. Local exponents, n, obtained for (tt^) using Eq. (4.8 1 



This procedure has the advantage that the unknown coefficient a is ehminated, but 
the disadvantage that it is sensitive to noise. Using the va lues of xq obtained from the 
regression method, the estimates of n obtained using Eq. (4.8), as applied to (m^), are 
shown in Figure [Sj It is clear that these estimates of n are sensitive to noise for x 
greater than, say, 6.5m. Never-the-less, reasonably constant values are found for the range 
2m < X < 6.5m, and averaging over this range provides a single estimate of the exponent 
n for each grid. Of course, this estimate varies slightly depending on whether data for 
(u^) or (g^) is analyzed, and on the precise range of x used to form the average. The 
range of exponents thus obtained for each grid is also shown in Table 1. It is reassuring 
that these values are close to those found by the regression method. 

4.3. Maximum decay range method 

In this method Xq is chosen to give the longest power-law range in x and then n is 
determined for that value of xq. The procedure works as follows. First Xmax is fixed at 
the value given by the regression method (at around 8m). Next, a series of different values 
of Xq are chosen, and for each xq the best fit value of n is obtained for a range of values 
of Xrnin- The valuc of Xq which gives the widest range of constant n between Xmin and 
Xmax is then deemed to be the correct one. The exponent n then follows. 

T he app lication of this method to the conventional grid is shown in Krogstad & David- 
(2010|) and demonstrated here by applying it to the (g^) data from msg2 in Figure 
the other methods, slightly different values of n are obtained for a given 
grid, depending on whether (it^) or (g^) are analyzed. The values obtained are given 
in ranges indicated in Table 1, where it can be seen that all three procedures generate 
similar exponents. 

There are three striking features of Table 1. First, all three grids yield decay exponents 
very close to the classical Saffman value of n = 6/5. Second, the conventional and multi- 
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Figure 9. Decay exponents, n, for the grid msg2 as a function of Xmin for different values of 

the virtual origin xq. 



scale grids produce almost identical results (to within experimental uncertainty). Third, 
these decay exponents are a long way from the n in the range 1.7 to 2.0 quoted in |Hurst| 



& Vassilicos (2007) for their four fractal cross grids. It is likely, therefore, that all three 
of our grids produce Saffman turbulence. 

We can check the hypotheses that we h ave Saffman turbulence if we pre-empt our 

shows (ul) £^ /U'^£q versus {x~'Xo)/£q for all 



discussion of length scales in § 5 Figure 



10 



three grids. Recall that, in Saffman turbulence, self-similarity of the large scales requires 
^2^3 = constant, as distinct from, say, Batchelor turbulence, in which u^£^ = constant. 
It is clear from Figure [T0| that there is an initial transient, which is more or less restricted 
to {x — xo)/£o < 100 {x < 2.7m), after which u'^£^ is indeed more or less constant in all 
three cases. The scatter in this data is largely a consequence of the difficulty associated 
with estimating £, as will be discussed in § [Sj 

Finally we consider the dimensionless dissipation coefficient A in Eq. (4.5), which 



is normally taken to be constant during the decay of isotropic turbulence. Assuming 
isotropy at the small scales, the viscous dissipation rate, e, can be written as 



3Au^ f du 

2— = i^'^\b-. 



(4.9) 



which allows us to estimate A from measurements of (^{dux/dx) u and £. The corre- 
sponding values of ^ are plotted in Figure[TT]for all three grids as a function of {x—xo)/£o- 
Again, there is some scatter, largely due to the difficulty in estimating £. Never-the-less, 
once the turbulence is fully developed, say for {x — xo)/£o > 100, there is a slow but 
steady decline in A which is consistent across all three grids. As noted earlier, this slow 




Multi-scale Grid Turbulence 



13 



8 
7 
6 
5 
4 
3 
2 
1 




o 
















• Conventional grid 
□ IVIulti-scale grid#1 
O IVIulti-scale grid#2 




; □ 

■ o - 








; 
















: o □ 
• o □ 


□ 












-<>-.0^<> 


• 

B • ---^-eP** 


• •• •• 






































100 200 300 



400 



500 



Figure 10. {ul)£^/U^£ versus {x — xo)/£o for all three grids. 

variation in A means that, even if we have SafFman turbulence, with u'^i^ = constant, 
we need not recover n = 6/5, as this value of n relies on A being strictly constant. 



5. The Streamwise Development of Length Scales 

We now turn to the length scales 77, A and £. The Kolmogorov microscale, define d as 
77 = (z^'^/e)^^'', can be determined from the isotropic estimate of dissipation in Eq. (4.9) 
or from the decay rate of (q^), both estimates giving virtually identical estimates for rj 
in this type of flow, as demonstrated by Krogstad & Davidson (20101. The streamwise 
development of the Kolmogorov and Taylor microscales, 77^ and A^, are shown in Figure 

m 



Noting that combining Eq. (4.5 I and (4.6 1 yields 

10i^(x- 



All nil 

we see that A^ should scale as A^ ~ {x — xq), which is indeed verified in Figure 
is, in effect, confirmation of a power-law form of energy decay. 

As noted in § [2| the integral scale, i, is defined in the usual way as 



12 



£ = 



{ux{x)ux{x + r 



)) 



dr = 



f{r)dr 



(5.1) 
This 

(5.2) 



where /(r) is the usual longitudinal correlation function and, in practice, {ux{x)ux{x + r)) 
is evaluated using Taylor's hypothesis. Ideally /(r) should decay monotonically to zero 
for large r, but in experiments /(r) almost always exhibits a weak oscillatory tail which 
persists for many multiples of £. This makes it difficult to evaluate f{r)dr accurately 
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Figure 11. Streamwise development of A obtained from Eq. (4.9 1. 
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Figure 12. Streamwise distributions of rj^ (filled symbols) and (open symbols). The 

straight line is ~ (s — xo). 




since it will depend on where the integral is terminated, and so the integral in Eq. (5.2) 
was only taken up to the first zero crossing, which introduces a small systematic error. 
As a result, there is some uncertainty in the calculated values of t. Our estimates of i 
for all three grids are shown in Figure [Ts} 

Evidently, £ exhibits considerably more scatter than 77 or A, but it is clear that the 
development of (. is similar for all three grids and follows a power law. Moreover, power- 
law energy decay with constant A demands t ^ ut, and so SafFman turbulence predicts 
£ ~ (a; — Xo)°'^, provided we ignore the slow decline in A. This is shown in Figure 13 for 
comparison and the data for all cases follow such a trend for {x — xo)/£ > 100. 



6. Spectral Development of the Turbulence 

We have seen that, as far as the behavior of the integral scales u and £ are concerned, 
the conventional and multi-scale grids behave in almost identical ways. However, it might 
be argued that, because energy is injected across the scales with a multi-scale grid, the 
evolution of the energy spectrum may exhibit non-classical features. So we conclude by 
comparing the spectra of the multi-scale and conventional grids. We shall focus on the 
usual one-dimensional spectrum, Fii{k), defined by the transform pair 

2 f°° 

-Fii(fc) = - / {ux{x)ua;{x + r)) cos{kr)dr , (6.1) 
71" Jo 

poo 

{ux{x)ux{x + r)) — I Fii{k)cos{kr)dk . (6.2) 
Jo 

Let us first consider the issue of the energy injection scales for the multi-scale grids. 



Figure 14 shows Fii{k) for these two grids, measured just downstream of the mesh (at 




X ~ I.Sto), and normalized by the Kolmogorov microscales, rj and v. The mesh sizes 
Mi and bar widths t^, are marked on the horizontal axis using the rule of thumb that 
k = tt/£. It is clear that in both cases the range of geometric scales spans most of the 
turbulent spectrum. 

Next there is the issue of whether or not the turbulence produced by the multi-scale 
grids behaves in the classical way, with the spectrum collapsing at high k when normalized 
by the microscales rj and u, and at low k when normalized by the integral scales £ and 



u. Figure 15 shows Fii{k) for grid msgl measured at four different streamwise locations 
and normalized by (a) and rj and v and (b) i and u. In case (a) there is excellent collapse 
of the spectra at high k and no collapse at low k, while in (b) there is reasonable collapse 
at low k and no collapse at high k. The same plots are given for rasg2 in Figure [T6| and 
again there is excellent collapse at high k when normalized by rj and u, and at low k when 
scaled by £ and u also in this case. Once again, we see that the turbulence produced by 
the multi-scale grids behaves in a classical manner. 

Finally we compare spectra at the same location generated by the three different grids. 
Figure 17 shows Fnik) for all three grids, normalized by the integral scales t and w, at 
(a) £ ^ {x — xq)/£q w 40 and (b) £ ^ {x — xo)/£o ~ 290. There is excellent collapse at 
both streamwise locations, confirming that the nature of the turbulence coming off the 
conventional and multi-scale grids is essentially the same. The same information may be 
examined in real space using the longitudinal structure function 

((Aw,)') = (^{u,{x + r) - u,{x)f) = 2 _ 2 {ujx)u,{x + r)) , (6.3) 

which, crudely speaking, a cts as a cumulative index of all energy held in structures 



smaller than scale r. Figure 



18 shows {^{Aux) y /2(u^) plotted against r/£ at the same 

two streamwise locations asTugure [TtI Again it can be seen that there is good collapse 
at both locations. 
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(a) (b) 
Figure 15. Fi\{k) for grid msgl measured at four streamwise locations and normalized by (a) 

rj and v, and (b) £ and u. 




7. Conclusions 

Our primary findings are two-fold. First, it seems that Saffman's decay law is reason- 
ably robust, since the energy decay exponents for all three grids are close to Saffman's 
classical prediction of n ~ 6/5. Second, the multi-scale grids used here produce almost 
identical results to the equivalent classical grid. In particular, all three flows exhibit 
remarkably similar streamwise distributions of Re\ (Figure [2| , flatness and skewness 



(figure 5), and dimensionless decay coefficient A (Figure 11 ). It is also worth noting that 
and the spectra for the multiscale grids exhibit classical Kolmogorov scaling, with E{k) 
collapsing on £ and u at low k, and on rj and v at high k. 

Our findings contradict those of some previous studies which report unusual behav- 
ior behind similar multi-scale grids, in particular, a very high energy decay exponent of 
around n ~ 2.0 and unusually high values of Re\. A decay exponent of n ~ 2.0 is par- 
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ticularly worrying as the theoretical maximum for n (assuming the dimensionless decay 
coefficient, A, is constant) is n = 10/7. However, these earher measurements were taken 
much closer to the grid where the flow exhibits initial grid-dependent inhomogeneities; 
inhomogeneities which, according to the present data, disappear further downstream. 



Appendix A. Classical decay exponents 



, including those of 


Kolmogorov 


Ossai & Lesieur 


2000 


Ishida et 



al. 2006 Davidson 2004 and Davidson 2009[ for more details.) We start with the result 
that, in fully-developed turbulence, E {k ^ Q) = Cmk"^, where Cm ~ u'^£™~^^ and Cm = 
constant for m < 4. Self-similarity of the large scales then demands u 



2fm+l _ 



constant 
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for TO < 4. When combmed with the zeroth law, Eq. |4.5[ this yields 

^ 1 + ^ , n^2 , 




n to m + o 



(Al) 



where the subscripts indicate values at < = 0. The smallest value of m is to = 2, since 
lower values require that the spectral tensor diverges as fc — > 0. On the other hand, the 
largest value is to = 4, as higher values spontaneous convert to m = 4 during an initial 
transient. Thus the classical range of decay exponents is, 6/5 < n < 10/7. 
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